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Abstract 

As an extension of our previous work on the holon pairing instability in the 
t-J Hamiltonian [Phys. Rev. B 66, 054427 (2002)], we examine the orbital 
symmetries of holon pairing order parameters in high Tc superconductivity 
by examining the energy poles of t-matrix. We find that both s- and d-wave 
symmetries in holon pair order parameter occur at low lying energy states 
corresponding to the higher energy poles of t-matrix while only the s-wave 
symmetry appears at the lowest energy pole and that this results in the d- 
wave symmetry in the Cooper pair order which is a composite of the holon 
pair of s-wave symmetry and the spinon pair of d-wave symmetry below Tc- 
Finally we demonstrate that there exists no time-reversal symmetry breaking 
in association with the Cooper pair order parameter. 
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I. INTRODUCTION 



At present the Cooper pair order parameter of dx^-y^-wave symmetry for high Tc cuprates 
is generally accepted [1]. However, recent tunneling measurements on YBa'^Cu^O'j^^ com- 
pound show that pairing symmetry changes from a pure (^3;2_2^2-wave symmetry to a mixture 
of d- and s- wave symmetries with the change of hole concentration [2-4] or with the varia- 
tion of applied magnetic field [2] . Cooper pair order parameters of complex mixing such as 
dx2_y2 +idxy and d^p^y^ -\-is break the time-reversal symmetry [2,4]. There have been some 
theoretical attempts to explain the observed thermal conductivity [5] in high Tc cuprates by 
resorting to the time-reversal symmetry breaking of the Cooper pair wave function [6,7]. On 
the other hand, the signatures for time-reversal symmetry breaking have also been found 
in the angle resolved photoemission spectroscopy [8] , neutron scattering [9] and /x-SR mea- 
surements [10] below the pseudogap temperature T*. Earher, Varma [11] proposed that the 
fourfold pattern of circulating current inside the Cu02 unit cell breaks the time- reversal 
symmetry. Chakravarty et al. [12] suggested that the time- reversal symmetry is broken by 
the d-density wave order which involves the circulating current around the Cu-0 bond. 

Recently we [14] reported the superconducting phase diagram by applying the Bethe- 
Salpeter equation to our earlier U(l) and SU(2) holon-pair boson theory [15] of t-J Hamil- 
tonian and obtained an arch shape structure of superconducting transition temperature in 
agreement with experiments. The da;2_y2-wave symmetry of Cooper pair can arise as a com- 
posite of the s-wave symmetry of the holon pair and the (i^2_j^2-wave symmetry of spinon 
pair. In the previous t-matrix study we briefly reported a discussion on the li-wave symme- 
try of the Cooper pair order parameter at and below T^. Here we present a detailed study 
on how the orbital symmetry of the holon pairing order parameters varies as a function of 
excitation energy below and above Tc- It is of great interest to examine whether the widely 
used t-J Hamiltonian can intrinsically predict the time reversal symmetry breaking by al- 
lowing a complex orbital mixing. The objectives of the present paper are two-fold. First we 
discuss the variation of the orbital symmetry for holon pair order parameter with excitation 
energy by examining the t-matrix pole. Second we examine whether there exists the time 
reversal symmetry breaking in association with the Cooper pair order parameter based on 
the U(l) and SU(2) slave-boson theories of the t-J Hamiltonian. 
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II. U(l) AND SU(2) SLAVE-BOSON REPRENSENTATIONS OF T-MATRIX 

BASED ON T-J HAMILTONIAN 



We write the t-J Hamiltonian, 

H^-t ^(ctcj> + C.C.) + J ^(Si • Sj - ^niUj), (1) 

where Sj • Sj — \ninj — — |(4j.4t ~ 4T^lj.)(^jT^iJ. ~ ^ji^i})- Here Sj is the electron spin 
operator at site i, Sj = ^cla'^apCip with a^p, the Pauh spin matrix element and rij, the 
electron number operator at site i, rii = cl^Ci^- In the slave-boson representation [16-21] 
the electron annihilation operator of spin a is written as a composite of spinon /g. (spin 
1/2 and charge object) and holon b (spin and +e object) operators. That is, — fa 
in the U(l) theory and Cq = ■^h^'^a in the SU(2) theory with a = 1, 2, where fa{i>) is the 

J2 



spinon (holon) annihilation operator in the U(l) theory, and ■^i = 
/ 



^ and -02 = ^ 
\ /2 / \ -fi 



and /i 



are respectively the doublets of spinon and holon annihilation operators 



V 



&2 

in the SU(2) theory. In the slave-boson representation, the Heisenberg interaction term is 
written Hj = -lY.hbjb]b\{fy]^ - fy]{){fj^h - Ink) the U(l) theory and Hj = 
-i E(l + /il/ii)(l + h]hj){fyi - flfj^Mih - fj2fii) for the SU(2) theory. In the above 
expressions the coupling between the charge (holon pair) and spin (spinon pair) degrees of 
freedom naturally arises from the composite nature of electron having both spin and charge. 

After decomposition of the Heisenberg interaction term and proper Hubbard Stratonovich 
transformations, the effective mean field Hamiltonian of holon is derived to be [15], for the 
U(l) theory (see Appendix A for a derivation), 

-^EHl'bWM-^Eblh, (2) 

where Xij =< fLfja + ^blbj > is the hopping order parameter and A{j =< fj^fn-fjji^ >, 
the spinon pairing order parameter, and for the SU(2) slave-boson theory, (see Appendix B 
for a derivation). 
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Ht-j,su{2) — ~ 9 X] hlUijhj + c.c. 



where Uij 



a,f3 



2 

E lAf/fcL&Jv^&.-^&ia - E h\K (3) 



/* 



is the order parameter matrix of hopping order, xa ^-nd spinon 

Ar* 

pairing order, A|j. with Xij =< fLfja + ^^(^li^ji - &]2^i2) > and Af^- =< /,i/,2 - /j2/a >• 

From the Bethe-Salpeter equation [14] , we obtain a matrix equation for the t- matrix, for 
the U(l) slave-boson theory (see Appendix C for a detailed derivation), 

E (<^k',k" - ^k',k"(q> Qo)) V',k(q> %) = - k), (4) 
k" 

where 

1 , ,,,, n(6(k")) + e/^-(--"+<^)n(6(-k" + q)) 
-.',."(q,^o) = ^-(k -k ) ,g„-(e(-k" + q) + e(k")) ' 

v{\s.) — — J|A/p7k, the momentum space representation of holon-holon interaction with 

7k = {coskx + cos ky) and n(e) = -^^t—i, the boson distribution function. Similarly, we 
obtain, for the SU(2) slave-boson theory (see Appendix D for a derivation), 

linn \ II II 

ct p ct p I \\ iCt p ap / 



E ( V,k"5a'a" - K'l" ^o) ) C'fk 9o) 

= ^(k' - k)<5„.„V/3> (6) 



where 

Si? ,(-k"+q) 



/ ^/ // ^11 



a /3 a /?' 



k' ,k" 



, , lv^„' ,",^(^a;(k)) + e n(E,,(-k +q)) 
(q'?o) = ^E^(k -k) r.^^.^j, X 



^4;^ ' zgo-(^«;(k") + ^^;(-k" + q)) 

t^a'a;(k")C/,',;(-k" + q)C/^^„(k")[/;,^„(-k" + q). (7) 

Here -Ei(k) — — /i and £^2(k) = — — A* are the energy dispersions of upper and lower 

bands of holons with £^k = t\/ (X7k)^ + (A/(/7k)^ and (/7k = (cos/ca; — cos ky). Uapiii) — 
( 

is the unitary matrix which diagonalizes the one-body holon Hamiltonian with 

yvi, Mk 



Mk = ~ and fk = ^^^^^^ + (Here sgn((/9fc) denotes the sign of (cosA;^ — 

cosky).) It is noted that the t-matrices, tk',k(q> ?o) and ^k'^"^(q' ^o) U(l) and SU(2) 



theories respectively are independent of the Matsubara frequencies /cq and ko, owing to the 
consideration of instantaneous holon-holon interaction, v(k' — k). 

The t-matrices are numerically evaluated from the use of Eqs.(4) and (6) at each tem- 
perature and hole doping concentration. The hopping and spinon pairing order parameters 
in Eqs.(2) and (3) are the saddle point values evaluated from the usual partition functions 
involving the functional integrals of slave-boson representation. [15] Using the matrix equa- 
tions (4) and (6), the poles of the t-matrices are searched for as a function of energy with 
q = 0, i.e., the zero momentum of the holon pair by using 

EVk'(?o,q = O)iy(go,k') = XW{qo,k), (8) 
k' 

for the U(l) theory, where W{qQ,k) is the eigenvector and A, the eigenvalue and 

E C'''''(5o,q=0)li^,y(go,k')=Aiy,,^(go,k), (9) 

k ,a ,/3 

for the SU(2) theory, where Wai3{qo: k) is the eigenvector associated with the SU(2) isospin 
channels a and p{a = 1,2 and P = 1,2 ) and A, the corresponding eigenvalue. The 
divergence of the eigenvalue at a given energy qo signifies the pole of the t-matrix at the 
energy. The value of pole go corresponds to the energy of holon pair corresponding to the 
eigenvector W{qo,\s.). 

For each eigenvector we compute the contribution of the s-, Px-, Py- and (i-r2_j^2-orbitals. 
The orbital symmetry is determined from numerical fitting to the computed result of the 

eigenvector, 

W{qo, k) = as(cos kx + cos ky) + Up^ sin kx + sin ky + ad{cos kx — cos ky) (10) 
for the U(l) theory and 

Wa/3{qo, k) = tts^afsicos kx + COS ky) + ttp^^afS sin kx + ttpy^afi sin ky + aa^afsicos kx - cos ky) (11) 

for the SU(2) theory. Here, ai represents the weight of the l-th partial wave of holon pair 
in the U(l) theory with orbital angular momentum oi I — s, Px, Py and d. For the SU(2) 
theory, ai^ai3 represents the weight of the l-th partial wave in the ba and 6/3 holon pairing 
channel with a,l3 = 1, 2. The poles and eigenvectors are obtained for three different hole 
concentrations : underdoped, optimally doped and overdoped cases with J/t = 0.3 and 
= 10 X 10 square lattice. 
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III. RESULTS 



In Fig. 1, we summarize our results for the energy poles and corresponding eigenvectors 
of the t-matrix. Above the on-set temperature Tc, there is no phase coherent holon pair 
bound state. As temperature is lowered to a temperature below from a temperature 
above Tc, the lowest energy pole changes from a positive to negative sign, that is, a phase 
coherent bound state with negative energy occurs [14]. The predicted bound state has the 
s-wave pairing symmetry. As shown in Fig. 1, the energies of the higher lying pairing states 
remains to be positive (unbound) even below Tc. 

In Table 1, we show the weight of each orbital for the eigenvector of the U(l) t-matrix 
at the lowest energy at temperature slightly above the on-set temperature {T = Tc + 1Q~H 
with the hoping integral). The eigenvector of the lowest energy corresponds to the stable 
channel of holon pairing below Tc [14]. As shown in Table 1, stable holon pairing occurs 
in the pure s-wave channel. The contributions of the higher partial wave (p- and d-wave) 
symmetries virtually vanishes; they are found to be less than 10~^^ times smaller than that of 
the s-wave contribution. The lowest lying holon pairing state remains to be the pure s-wave 
(having no mixing with other orbitals) with decreasing temperature from a temperature 
above Tg. It is now clear from this prediction that the Cooper pairing occurs with the pure 
dx2-y2 symmetry resulting from the composite nature of the d^'^-yi-wave symmetry of spinon 
pairing and the s-wave symmetry of the holon pairing [22] . 

With the SU(2) t-matrix, there exist two eigenvectors corresponding to two bosons bi and 
62- The weight of each orbital for the first eigenvector is shown in Table 2 and the second, 
in Table 3. As shown in Table 2, the first eigenvector is of pure s-wave pairing symmetry 
with no phase difference between the bi-bi and 62-^2 pair scattering channel. As shown in 
Table 3, the second eigenvector shows a dominance of the s-wave pairing symmetry with the 
phase difference of tt between the 61-61 and 62-62 channels and a negligibly small contribution 
from the (i-wave symmetry in the 61-62 channel. The weights of other orbital contributions 
are negligibly small as shown in the Tables. It can be readily checked that only the second 
eigenvector in Table 3 leads to non-vanishing Cooper pair order parameter, while the first 
one in Table 2 results in vanishing Cooper pair order parameter [23]. The SU(2) theory 
differs from the U(l) theory in that the phase fluctuations of order parameters are taken 
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into account by allowing the 62-boson [21]. There are two eigenvectors of t-matrix associated 
with the same energy pole in SU(2). However, there is only one stable s-wave channel of 
holon pairing with no other orbital mixing which leads to the non-vanishing Cooper pair 
order parameter for both the U(l) and SU(2) theories. It is now clear that the composition 
of the s-wave channel of the holon pair and the d-wave channel of the spinon pair leads to 
the pure d^2_2^2-wave symmetry of the Cooper pair order parameter [1]. Thus there exists 
no time-reversal symmetry breaking allowing no complex orbital mixing. 

For the higher energy poles, the s-, px-, Py- and (i^2_j^2-wave states are predicted to occur 
as shown in Tables 4 and 5. Interestingly, the eigenvectors occur in a sequence of certain 
pattern as the energy pole of the t-matrix increases both above and below Tc as shown in 
Fig. 1 : s-wave state occurs in the lowest energy; the next four high lying energy states are 
of the S-, Px-, Py- and dx2^y2-waye symmetries with nearly the same energy (that is, near 
degenerate); the next higher lying energy states are of s-, px- and Py-wave symmetries with 
nearly the same energy and so on. Below we will discuss how the above sequential pattern 
arises in association with momenta available in the intermediate scattering states of holon 
pair . 

First we analyze the locations of the poles. The poles of the t-matrix occur at discrete 
energies for both bound and unbound states of the holon pairs for finite lattice. With 
increasing lattice size {L — > 00), the unbound states become continuum as the allowed 
momenta become continuum. The energy of holon pair with momenta (k and — k) is the 
sum of kinetic energy and interaction energy. The magnitude of the interaction energy 
{V ~ J|A/p) is relatively small as compared to that of the kinetic energy {K ~ tx) with 
the ratio ofV/K < 0.1. Thus the kinetic energy dominant ly determines the energy of holon 
pair. Indeed, the predicted locations of poles in Tables 4 and 5 are close to the kinetic energy 
of holon pair with momenta k = (0, 0), (27r/L, 0), (27r/L, 27r/L) and (27r/L, 0) [24]. 

Here we pay attention to the intermediate scattering states of holon pair. This will allow 
us to examine how an initial holon pair state with an l-th partial wave states {I — s,Px,Py 
or d) can be scattered into intermediate states having momenta ki and — ki. The transition 
amplitude from the initial state to the intermediate state is, to first order, 

Ti^(ki,-ki} = ^v(ki - k)wi{k) 

k 
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= - J| A/p [cos(A;ia; - k^) + cos{kiy - ky)] ^^(k), (12) 

k 

where wi{k) is the eigenvector with an l-th partial wave (orbital). The Feynmann diagram 
for this process is displayed in Fig. 2. 

For the initial s-wave pairing state ^^(k) = {cosk^ + cos ky), the transition amplitude is 
obtained to be 

7;^{ki -ki} = --^1 A/P J2 [cos(/cia; - k-^) + cos{kiy - ky)] (cos k^^ + cos ky) 

k 

= — ^J|A/|^(cos /ci^ + COS /ciy). (13) 

This amplitude does not vanish as long as the momentum of the intermediate state (ki) 
satisfies the condition of (cos /ci-c + cos k\^ ^ 0. In other words, the s-wave state has a 
non-vanishing transition amplitude to the intermediate state with momenta ki and — ki if 
(cos kix + cos kiy) 7^ 0. Thus, the s-wave scattering channel allows poles at the intermediate 
state energies of holon pairs with momenta ki and — ki as long as (cos ki^ + cos kiy) ^ 0. For 
example, the s-wave eigenvector occurs at poles near 2e(0, 0), 2e(27r/L, 0), 2e(27r/L, 27r/L) 
and 2e(47r/L, 0) as shown in Table 4. This is because (cos/cix + coskiy) does not vanish 
for momenta, (0,0), (27r/L, 0), (27r/L, 27r/L) and (47r/L, 0). Similarly, one can readily un- 
derstand why p- and c?-wave scattering states (channels) occur only at certain energy levels 
[25]. For higher angular momentum channels such as Wi{k) = cos2kx or wi{k) = sin2A;a;, 
the transition amplitude in Eq. (13) identically vanishes. Therefore, t-matrix can not allow 
poles in the scattering channels corresponding to these higher orbital angular momentum 
states; only the s, Px, Py and dx2-y2 state can occur. 

The above analysis with the U(l) theory is readily apphed to the SU(2) case to find the 
location of poles and the corresponding eigenvectors. The only difference is that in SU(2) 
the number of eigenvectors for each pole is doubled owing to the presence of two degenerate 
bosons bi and 62- These two eigenvectors are found to be degenerate. One is the eigenvector 
with no phase difference between the bi — bi and 62 — ^'2 pairing channel. The other is the 
eigenvector with the phase difference of tt. Interestingly, there is a small component of 61 — 62 
boson pairing for the case of the latter. The weight of this off-diagonal pairing component 
is negligibly small in the lowest energy pole (e.g. ~ 10""^ for u = 0.0097 with x = 0.07), 
but become appreciably large at higher energy poles (e.g. ~ 0.2 for uj/t — 1.1236 with 
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X — 0.07). On the other hand, the weight of the bi — 62 boson pairing state decreases with 
the increase of hole concentration as shown in Table 5. This is due to the decreasing trend 
of the spinon pairing order parameter Af with increasing hole concentration. It is of note 
that the off-diagonal element between bi and 62 bosons (or vice versa) in the Hamiltonian 
is proportional to the spinon pairing order parameter in Eq.(3). 

It is also noted that there occurs not only s- and d- wave states, but also p- wave 
state for holon pairing. The p- wave state is not allowed because it has the odd value of the 
orbital angular momentum and violate the symmetric waves. With the inclusion of exchange 
channel only s- and d- wave states appear in association with eigenvectors in the t-matrix 
as shown in Fig. 3. 

IV. SUMMARY 

In the present study, we investigated the orbital contributions to both U(l) and SU(2) 
t-matrices both above and below Tc and examined the lowest energy and higher energy 
poles respectively by applying the Bcthe-Salpeter equation to the U(l) and SU(2) holon- 
pair boson theory. We showed that at the lowest energy (below Tc) there is only one holon 
pairing bound state which gives rise to the non- vanishing Cooper pair order parameter for 
both the U(l) and SU(2) theories. We found that the holon pairing at and below Tc is 
made of the pure s-wave with no complex orbital mixing and thus the Cooper pair of d- 
wave symmetry as a composite of the s-wave holon pair and the d-wave spinon pair. As a 
consequence there exits no time reversal symmetry breaking. Thus we argue that symmetry 
breaking is caused by a different origin since the complex orbital mixing such as dx2^y2+idxy 
and dx2^y2 + is did not appear. 
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APPENDIX A: DERIVATION OF THE U(l) HOLON HAMILTONIAN 



The t-J Hamiltonian of interest is given by, 

H ^-t {cl^Cj^ + c.c.) + J ^ (Si • Sj - jUirij) 

<i,j> <i,j> 

-I^J^claCia (Al) 
i 

and the Heisenberg interaction term is rewritten 
Hj^J ^ (Si • S^- - jriiTij) 

<i,j> 

= - 9 E (4ic]t - c\Ai){c3^Cii - CjiCit). (A2) 

^ <i,3> 

Here t is the hopping energy and Sj, the electron spin operator at site i, Sj = \c^ia(^ apcip with 
(To,/3, the Pauli spin matrix element, rii is the electron number operator at site i, rii — cl^Cia- 
jj, is the chemical potential. 

In the U(l) slave-boson representation [17-20], with single occupancy constraint at site 
i the electron annihilation operator Qo- is taken as a composite operator of the spinon (neu- 
trally charged fermion) annihilation operator /jg. and the holon (positively charged boson) 
creation operator 6|, and thus, Cjo- = fiabl- Rigorously speaking, it should be noted that 
the expression Qo- = fef/io- is not precise since these operators belong to different Hilbert 
spaces and thus the equality sign here should be taken only as a symbol for mapping. Us- 
ing — fiabl and introducing the Lagrange multiplier term (the last term in Eq.(A3)) to 
enforce single occupancy constraint, the t-J Hamiltonian is rewritten, 

H^-tY. {{flk){blfja) + c.c.) 

<i,j> 

E b^b,b]bi{fy], - lyimh - ink) 

<ij> 

-1^ E flhhabl - I E Hb\h + flh. - 1). (A3) 

i i 

The coupling between physical quantities A and B is decomposed into terms involving 

fluctuations of A, i.e., {A— < A >) and B, i.e., {B— < B >), separately uncorrelated mean 
field contribution of < A > and < B > and correlation between fluctuations of A and B, 
that is, {A- < A>) and {B- < B >) respectively; AB ^< B > A+ < A > B- < A >< 
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B > +{A- < A >){B- < B >). Setting A = hihjh]h\ for charge (holon) contribution 
and B = [flj]^ - fy]i){fj^fii - fjifi]) for spin (spinon) contribution, the Heisenberg 
coupUng term can be decomposed into terms involving coupling between the charge and 
spin fluctuations separately, the mean field contributions and coupling (correlation) between 
fluctuations (charge and spin fluctuations). Using such decomposition of the Heisenberg 
interaction term for Eq.(A3), we write the U(l) Hamiltonian 

E(/U4&^ + c-c.) 

<i,3> 

-i E [((44-44)(/.T/a-/.i/.T)M%' 

^ <i,j> 

-(b.b,b\bl){ifyh - fUnU^h - hik)) 

+{kbjb]bl - (bibjblbl)) 

-1^ E + b\bi) - I E Wlfia + b% - 1). (A4) 

i i 

Noting that [6^, 6^] = 5ij for boson, the intersite charge (holon) interaction term (the 
second term) in Eq.(A4) is rewritten, 

-^((44t - fyiMifil - fM)bibjb]bl 

= < |Af,.| >2 (l + blh + b]bj + blb]bjb,) , (A5) 

with Af — fjifii — fjifih the spinon pairing held. The third term in Eq.(A4) represents 
the intersite spin (spinon) interaction and is rewritten, 

-^{bibM){fyj, - fyimfii - fM 

- -^(44 - fyiMifri - fjM^ (A6) 

where Jp — J(l+ < > + < hjbj > + < blbibjbj >) or Jp — J(l — x)^ with x, the uniform 
hole doping concentration. The fourth term in Eq. (A4) is written, 

^{bibjb]bi){{fyj^ - fyimfii - fM) 

= - (l+ < bjbi > + < b]bj > + < bjbiblbj >) < |Af/ > . 

(A7) 
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The intersite spinon interaction term in Eq.(A6) is decomposed into the direct, exchange 
and pairing channels [20], 



J n 



'Wifilfjl ~ filfjl){fjlfil ~ fjlfil) 



^ k=l 
^ Vd + Ve + Vp 



(A8) 



with (7° = /, the identity matrix and cr^'^'^, the Pauh spin matrices, where vd, ve and vp 
are the spinon interaction terms of the direct, exchange and pairing channels respectively 



[15], 



° k=0 



Jp 

T 

J r, 



Here (j° is the unit matrix and cr^'^'^, the Pauh spin matrices. 
Combining Eq.(A5) and Eq.(A7), we have 
J 



(A9) 

(AlO) 
(All) 



< \A{/ > (l + blk + h]h, + h\h]hjh^ 



+^ < |Af/ > (l+ < blbi > + < b]bj > + < b%b]bj >) 

K? > blb]b,k + ^ < |Af/ X blkb]b, > 



- < blk >) + (b]bj- < b]bj >) 



(A12) 



Collecting the decomposed terms Eq.(A5) through Eq.(A7) in association with Eqs.(A8) 
through Eq.(Al2), we write 

<i,j> 



<ij> 



4 



k=0 



12 



J 



<id> 



J 



-2 E W 

^ <i,j> 



f |2 



[blbi- < h\hi >) + (h]h,- < b]bj > 



(A13) 



where we considered < |A|^ p >= |A{^ p and ignored the fifth term in Eq.(A4). 

As arc shown in Eqs.(A9) through (All) the spinon interaction term is decomposed 
into the direct, exchange and pairing channels respectively. Proper Hubbard- Stratonovich 
transformations corresponding to these channels and saddle point approximation leads to 
the effective Hamiltonian, 



^eff = X E 



4t 



<i,j> 



J 



77 E Hm^jbi 



<hi> 



+f EKr-Aj(4iJT-44)-- 

^ <i,j> 

+|^EE[(p5-r-p;(//-'/0 

^ <i,j> 1=0 

At" 



+■ 



J, 



E ip\b,){b]b,) 



P <i,j> 



J 



+- y: \Ai\' < bib, X b]bj > 



J 



- E |Af/ [{b\h- < blk >) + {bib,- < bib, >) 



<^,J> 



-/^ E + 4^.) - ^ E HfU- + ^I^^ - 1)' 



(AM) 



where A^^. =< 6^6,- >, Xij =< flfja + f/i^j >, Af^- =< /,-^/^ - fjj^^ > and p,^ =< 
\fl(T^ fi > are proper saddle points. 

We note that = | < f^a^fi >=< -S^ >= for Z = 1, 2, 3, = i < /t/i^ >= i for Z = 
for the contribution of the direct spinon interaction term (the fourth term) . The expression 
{blbi){b\bj) in the fifth term of Eq.(Al4) represents the exchange interaction channel. The 
exchange channel will be ignored owing to a large cost in energy, U ^ ^ [20] [21]. The 
resulting effective Hamiltonian is 



(A15) 
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where H^'^ represents the the saddle point energy involved with the spinon pairing order 
parameter A-'^ and the hopping order parameter x, 



(A 16) 



H'^ is the holon Hamiltonian, 



<i,j> 

J 



<i,3> 



i 

where /x^ — %\ + ^ J2j=i±x,i±y I^^ P ^ spinon Hamiltonian, 



(A17) 



<i,j> 



-y E [A{;{f,,M-f,M + c.c. 



<^,3> 



E/^f(/l./^--(i-^))' 



(A18) 



where fj,{ — //(I — x) + iXi. 



APPENDIX B: DERIVATION OF THE SU(2) HOLON HAMILTONIAN 

The SU(2) slave-boson representation of the above t-J Hamiltonian reads [15] 

<i,j>a 

+ (/2i/lj - flif2j){b\jb2i + b\ib2j) + c.c. 

j:i^ + hlhi)il + h]hj)x 

^ <i,j> 

i 

i 

+i>^f\flihi - f2iA + blbu - 



(Bl) 



where A 



(1),(2),(3) 



are the real Lagrangian multipliers to enforce the local single occupancy con- 



straint in the SU(2) slave-boson representation [21]. Taking decomposition of the Heisenberg 
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interaction term above into terms involving charge and spin fluctuations separately, uncorre- 
lated mean field contributions and correlated fluctuations, i.e., correlations between charge 
and spin fluctuations as in the U(l) case, the SU(2) Hamiltonian is rewritten, 

<i,j> 

+{fi2fji - fiifj2){blibi2 + blibj2) + c.c. 

E [((/i^/ji - fy]2){foih2 - /,2/a))(l + h\hi){l + h]h,) 

^ <i,j> 

+ ((1 + hlh,){l + h%)){fyi - fy],){f^,h2 - fj2fil) 

-((1 + hlh.){i + h%)){{fyi - Afj^Mif^^ - fM) 

+ ((1 + hlhi){l + h]hj) - ((1 + hlhi){l + h]hj))) X 

{ifkfjl ~ filfj2)ifjlfi2 — fj2fil) — ((/e/jl ~ filfj2)ifjlfi2 " /j2/il))) 
i 

- EO^fVii/^L + bih^) + i\f\fi2fn + bihn) 

i 

+iAr - fi2f!2 + blbn - bib,,)) . (B2) 
The intersite holon interaction (the third term in Eq.(B2)) is rewritten, 

-^{(fyji - /A4)(/.i/^2 - /,2/a))(l + h^^^){l + h]h,) 

I A// > (1 + hlh, + h]hj + hlh,h%) 
~ + blk^ + b],b,o. + blb]f,b,0k^) 

(B3) 

where A^- = (^fjifi2 — /j2/ii) is the spinon singlet pairing order parameter. The intersite 
spinon interaction (the fourth term in Eq.(B2)) is rewritten in terms of decomposed Hatree- 
Fock-Bogohubov channels in the same way as in the U(l) case, 

-y (/^L/Ji - fkfj2)ifnfi2 - fj2fii) ^vn + VE + VP, (B4) 

where vd, ve and vp are the interactions corresponding to the direct, exchange and pairing 
channels respectively. 
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J, 



1=0 



Ve = -^{{flfja){flfia) -Ui), 



J, 



Vp = 



^{fhfjl ~ filfj2){fjlfi2 — fj2fil)- 



Ft 



(B5) 

(B6) 
(B7) 



Here cr" is the unit matrix and cr^'^'^, the Pauh spin matrices. The fifth term in Eq.(B2) is 
written, 

((/i2/jl ~ filfj2)(fjlfi2 — fj2fil 

[{1 + h^,h,)){{l + h%)) X 



ifLfjl - filfj2)){{fjlfi2 - fj2fil 

^|AJ|2(1+ < hlhi > + < h]hj > + < h\hi >< h]hj >). 



ft t\ 



(B8) 



We introduced ((/i/ji - /A4)(/,i/.2 - /,2/.i)) ^ {(/^/ji - /A4))((/.i/.2 - /,2/.i); 
|A{/ and ((1 + h%){\ + ^ ((1 + h\K))l^{\ + 

By introducing the Hubbard-Stratonovich fields, , Xij and Ajj for the spinon direct, 
exchange and pairing order shown in Eqs.(B5), (B6), (B7), we rewrite the effective Hamil- 
tonian for Eq.(B2), 



2t 



H¥u\2) [lX..r - X^AfUj + t(4&.1 - ^'l2^'2.)} - C.C 

Or) 



J 



<i,]> 

/|2kt At 



E blh]pho.hp 



<i,j> 



<i,j> P 
<i,j> 1=0 

+ y E [(^Ii^.i - b]2h2Wnba - h\M 



P <i,j> 



^ E H\'[hlh, + h\h~ < h]h, >-< h\h, >] + ^ E K 

<i,j> <i,j> 



/ |2 2 
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i i 



(B9) 



where Xij =< flfja + f - &l2&i2) >, A^^- = ((/ii/i2 - fi2fji) - -j;ip^i2bn + b]2hi, 

Mj - -j{hx)^%;i2^ ^i^'^ ^ii;i2 = (^L^ji + b]2hi) and a;, the hole doping rate. 

To simplify the effective Hamiltonian, we rearrange each term in Eq.(B9). For the 
paramagnetic state, we obtain p\ — \{flcr^fi) —< S\ >— for Z = 1,2,3 and p° = | < 
ficficr >— I for Z = 0. The one-body holon term (the sixth term) in the above Hamiltonian 
is incorporated into the effective chemical potential term. By setting Aij = A{j — jj^^z^ Xij -,12 
where A{j =< /ii/j2 — fi2fji > and xtj-u = (^i2^ii + ^j2^«i)) rearrange the third term 
and the second term in the bracket of the fifth term to obtain the effective Hamiltonian, 



<i,3> 



<i,j> 



— c.c. 



<«J> 



-f E te(/U.) + c.c. 

<i,j> 

I E H\%bhb,A, 

^ <i,j> 

-t ^ K*(/.i/^2-Wa) + c.a 

^ <i,j> 
i 

-EN(/i^4 + 44)+^A?(/^2/^i 



bi2bii) 



^ E (^f. - (/.i/^2 - Mn)) X^;;i2 - c.c. 



<^,3> 



2J, 



E |4-;i2 - iblbji + 6]2&ii) 



P <i,j> 



+-r E (^Ii^Ji - b]2bi2)ib],ba - b^bj^), 

"^P <i,j> 



(BIO) 



where /^^ 



/ 12 
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We neglect correlations between the fluctuations of order parameters. This is because 
correlations between the spin (spinon pair) and charge (holon pair) fluctuations are expected 
to be small as compared to the saddle point contribution of the order parameters (the flrst 
and second terms in the above Hamiltonian) particularly near the pseudogap temperature 
and the bose condensation temperature. However, individual fluctuations of the spinon pair- 
ing and holon pairing order parameters are not ignored. We also neglect the fluctuations 
of order parameters (ninth and tenth terms) in Eq.(BlO). The ninth term represents the 
fluctuation of spinon pairing order parameter and we neglect it owing to its vanishment as 
an expectation value. The tenth term represents the fluctuations of holon hopping order 
parameter and is negligible at low temperature as its fluctuations die out in the low temper- 
ature regions where pairing order parameters A-'^ begins to open. Owing to the high energy 
cost involved with the Coulomb repulsion energy the exchange interaction terms (the last 
positive energy terms) will be ignored [20]- [21]. We then obtain the mean fleld Hamiltonian, 
H = H^'^ + H'^ + , where H^'^ is the saddle point contribution of order parameters, x 
and A^, 

^5%) = ^ E lAjiv + :| E [I Air + \\x^^' + 1], (Bii) 

^ <i,j> ^ <i,j> ^ ^ 

is the holon Hamiltonian, 



<i,j> 



-A{^{b],k, + blbj2)] -c.c. 



<^,3> 



-Y:[l^'^ihlhi -x)+ iXliblbl) + i>Hibi2ba) 

i 

+iX',{blbn-bl2bi2)], (B12) 



and , the spinon Hamiltonian, 

4 

2" 



Hiu(2)--i E [xUfUj.) + C.C. 

<i,j> 



-y E HifjiU-fj2fii) + c.c: 

<ij> 



(B13) 
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where Xij =< flfj^r + f - &5-2&i2) >, =< fjifi2 - fj2fii > and fx^^ ^ -^^ - 

2 ^j=i±x,i±y Taking the saddle point value for the Lagrangian multiplier fields A' = 

[15], we obtain the holon Hamiltonian, Eq.(3). 



APPENDIX C: DERIVATION OF EQS.(4) AND (5) 

Introducing a uniform hopping order parameter, Xji — Xi and a o?-wave spinon pairing 
order parameters, = ±Ay(the sign +(— ) is for the ij link parallel to x (y)), we obtain 
the energy-momentum space representation of the action for the U(l) theory, 

k,a;„ 

J- ^ ^ 7;(k-k>(-k' + q,-a;;, + i/„)t&(k',u;;,)t6(k,u;„)6(-k + q,-a;„ + i/„), (CI) 



where 6(k, a;^) = ^= /o^ c^tZ^i e*^'^"'^~^''"'^6j(r), ^(k) = — 2tx7k — the energy dispersion of 



holon with 7k = (cos fc^: + cosky) and f (k' — k) = — J|A/p7k, the momentum space repre- 
sentation of holon-holon interaction. N is the total number of sites for the two dimensional 
lattice of interest and /3 = -^^j the inverse temperature. Considering the ladder diagrams 
for the holon-holon scattering, we obtain the Bethe-Salpeter equation for the t-matrix, 

< k' , —k' + q\t\k, —k + q >u{i)— v{\<i — k) 



7^ Yl '"0^ - ^')9{k")g{-k" + q) 



II ,11 



X 



<k ,-k +q\t\k,-k + q>u{i), (C2) 

where g{k) —< 6(k, a;„)6(k, >, the holon Matsubara Green's function and k — 
{cun, kx, ky) is the three-component vector of the energy and momentum. Using the fact 
that the holon-holon interaction v{\<i — k") is frequency independent, we sum the Matsub- 
ara frequency /cq. 



pY.9{k)g{ k +5)-^E,^;^'_^^k")z(-^o+go)-6(-k" + q) 

n(e(k")) + eP<-^'+M<-^' + q)) ^^3) 



igo-(e(-k" + q) + e(k")) 
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Inserting Eq.(C3) into Eq.(C2) and defining t^' ^(q, qo) =< k\ —k' + q\t\k, —k + q >u{i): we 
obtain the Bethe-Salpeter equation for the t-matrix, 



J2 ((^k',k" - "ik',k"(q, Qo)) ^k",k(q, Qo) = v{]^ - k), 



(C4) 



where 



N(3 



, - 

y iko - e(k") i{-ko + qo) - e(-k" + q) 



1_ , » n(e(k")) + e^^^-'^ +q)n(e(-k" + q)) 
n""^ ^ ^go-(e(-k" + q) + e(k")) 



(C5) 



with n(e) = ^3733;, the boson distribution function. 

APPENDIX D: DERIVATION OF EQS.(6) AND (7) 

Introducing a uniform hopping order parameter, Xji — X ^ d-wave spinon pairing 
order parameters, AJ^ = ±Aj(the sign +(— ) is for the ij hnk parallel to x (y)), we obtain 
the energy-momentum space representation of the action for the SU(2) theory, 

/ 



Ssu{2) = H (^l(k,^^n)^&2(k,u;„)^) 



-iujn + tX7k - A* 



61 (k, Un) 

62 (k, UJn) 



+^ ^(k-k')6^(-k' + q, -u;^ + i/„)^6„(k',u;^)^6c«(k,u;„)6^(-k + q, + 

k,k',q,a,/3a)„,a;^,fc'„ 

(Dl) 



where 6Q,(k, cj„) 



1 f,« 



So d^Y.i^^^'^^'^ ^'^^^haij) and (/9k = (cosA;^; — cos ky) with «, 



1, 2), two components of SU(2) holon. Considering the ladder diagrams for the holon-holon 
scattering, we obtain the Bethe-Salpeter equation for the t-matrix. 



< k\a';-k' + q,p'\t\k,a;-k + q,/3 >su(2)^ v(\i - ^)5^' ^Jp' 

J7S E ^(k' - ^')9a'a"{^')9p'0'{-^' + Q) 



k ,a ,p 



X <k ,a ]—k +q,P \t\k, a; -k + q, (3 >su{2), 
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(D2) 



where Qafiik) —< 6Q(k, a;„)6;3(k, a;^)^ >, the holon Matsubara Green's function. 

To calculate the holon Matsubara Green's function in the diagonalized basis of the one- 
body term in the action (Eq.(Dl)), we introduce the unitary transformation for the holon 
field, 



^ 62(k,a;„) J 



U{k) 



b'2{k,UJn) J 



(D3) 



where 



Uafsik) = 




(D4) 



with Mk = - Vi, 



V2 



1 + ^f^, Eik) = t^ixi^y + {Afcp^y and b'^iKuj^), 



the quasi-holon field. Then the one-body action in Eq.(Dl) becomes 

SsUi2) = 2^ (bi{k,UJny,b2[k,UJny) 



(^'l(k,'^n)^&2(k,^^r^)^) 



-iuJn + Ei{k) 




^ ^b[{k,u;n)^ 



-iu;n + E2{k) J 



\02 



(D5) 



where the quasi-holon energy is -E'i(k) = (X7k)^ + (A/(/7k)^ — A* E2{k) — 

— iy''(X7k)^ + (Afif-kY — jji. Prom the action in Eq.(D5) we readily obtain the holon Mat- 
subara Green's function, 



Qapik) =< ba{k,Un)bl{k,UJn) > 

= U^a'i^) < b'^,{k,u;„)b){k,u;„) > (t/(k)t 







iw„+E2(k) ) ^0 



{U{k)%'^, 



(D6) 



where -E'i(k) and -E'2(k) are the holon quasiparticle energy in the upper and lower band 
respectively and sgn{(fik) denotes the sign of (cos/c^ — cos ky). It is of note that there are 
two energy bands owing to the two kinds of holons in the SU(2) theory. 

Using the fact that the holon-holon interaction v{k! — k") is frequency independent, we 
sum the Matsubara frequency /cq in Eq.(D2), 
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^Y.9a'a"ik")90'(i"{-k" + q) 



n{E^,{k")) + e/^^^i(-'-"+-^)n(ig,,(-k" + q)) 

zgo-(%(-k" + q)+E«,(k")) ■ ^ ^ 

Inserting Eq.(D7) into Eq.(D2) and defining k'^^^*^' ~'~ ^' /^'l^l^' + 

9,/^ >sf/(2), we obtain the Bethe-Salpeter equation for the t- matrix, 



I I II Ji 



k ,a 

= ^(k' - k)5„,„5^,^, (D8) 

where 

t^a'a;(k'>;,'^;(-k" + q)C/^^„(k")C/t,^„(-k" + q). (D9) 
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with i and j are the nearest neighbor sites. Here, wc used the SU(2) representation [21] 
of electron annihilation operator, that is, Cj| = -^{blifn + blf},) for spin up electron and 
— -^{b\ifi2 — blfii) for spin down electron. In the real space, the first eigenvector in 
Table 2 is written < bl^b]^ >— (5a,i5/3,i -l- (^0,2(^/9,2) iorj — i±xorj — i±y. Combining 
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in a non- vanishing Cooper pair order parameter with pure dx2_y2 symmetry. The first 
two terms in Eq.(DlO) yield the d-wave symmetry of Cooper pairs as a composite of the 
s-wave symmetry of holon pair (< biabja >) and the d-wecve symmetry of spinon pair 
(< crfiafj-a >)■ The last two terms in Eq.(DlO) also lead to the d-wme symmetry of 
Cooper pair owing to the d-wave symmetry of the 61-62 holon pair (< biibj2 >) and the 
uniform phase of the spinon hopping order parameter (< fl^fja >)■ 

[24] For instance, the lowest energy pole (O.Olllt) for x = 0.04 in the U(l) theory ap- 
proximately equals the kinetic energy of holon pair with momenta k = (0,0), i.e., 
2e(0, 0) = O.Oldt with t, the hopping energy as is shown in Table. 4. The energies of 
the next four higher lying states are predicted to be reasonably close to the kinetic 
energies of holon pair corresponding to momenta k = (±27r/L,0) or (0, ±27r/L) (i.e., 
e(27r/L,0) + e(-27r/L,0) = e(0,27r/L) + e(0, -27r/L) = 0.724t) as shown in Table 4. 
(Here 27r/L is the smallest possible nonzero momentum of or ky for L x L square 
lattice.) Similarly the energies of the next three higher lying states are also close to 
the kinetic energy, e(27r/L, 27r/L) -|- e(— 27r/L, — 27r/L) = 1.432i; the next four higher 
lying states have energies which are approximately the same as the kinetic energy, 
e(47r/L, 0) + e(-47r/L, 0) = e(0, Att/L) + e(0, -Att/L) = 2.578t. Similar arguments hold 
true for other hole concentrations in Table 4. The SU(2) results are shown in Table 5, 
reveahng a trend similar to the U(l) case above. 

[25] For the p^-wave pairing state, wiCk) = sinkx, the transition amplitude is similarly 
obtained to be 

7;,^{ki,-ki} = -yJ|A/|^sin/ci^. (Dll) 

For sin kix — 0, the transition does not occur from the pairing state of Px orbital to 
any intermediate state with momenta ki and — ki. This explains why there is no eigen- 
vector of Px or py at the lowest energy pole which comes from the holon of momentum 
(0,0). Similarly, for the d-wave pairing state, wi{k) = {coskx — cos ky), the lowest order 
transition amplitude is given by, 

Trf^lki -ki} = -^J\Af\^{coskix - coskiy). (D12) 
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The transition amplitude vanishes if kx — ±ky. This imphes that there is no d-wave 
eigenvector at the pole near the energy of 2e(0, 0), 2e(27r/L, 27r/L). This is because the 
rf-wave pairing state can not make transitions into intermediate states having momenta 
ki and -ki where ki = (0,0) or (27r/L, 27r/L). 
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TABLE CAPTIONS 



Table 1 Weights of the 

Px't Py~ ^-iid ci-wave contributions to the holon pairing order 
parameter corresponding to the lowest energy in the U(l) holon-pair boson theory. 

Table 2 Weights of the s-, p^-, py- and rf-wavc contributions in each ba-bp scattering channel 
with a, /3 = 1, 2 to the first holon pair order parameter corresponding to the lowest 
energy in the SU(2) holon-pair boson theory. 

Table 3 Weights of the s-, p^-, Py- and c^-wave contributions in each ba-bp scattering channel 
with a, j3 = 1, 2 to the second holon pair order parameter corrsponding to the lowest 
energy in the SU(2) holon-pair boson theory. 

Table 4 The orbital state and energy at various hole doping x and temperatures {T/t — 
Tc/t + 0.001) for underdoping {x = 0.04, Tc/t = 0.034), optimal doping(a; = 0.07, 
Tc/t = 0.044) and overdoping [x = 0.1, T^/t = 0.041) respectively with the U(l) 
holon-pair boson theory. The calculations are done on the = 10 x 10 lattice with 
the use olJ/t — 0.3. Here, s denotes (cos A;a;-|-cos /Cj,); d, {coskx — cosky); Px, sinkx 
and Py, sin ky. 

Table 5 The orbital and energy at temperatures T /t = T^t + 0.001 for underdoping [x — 
0.07, Tjt = 0.027), optimal doping(a: = 0.13, Tjt = 0.037) and overdoping 
{x — 0.19, Tc/t — 0.031) respectively in the SU(2) holon-pair boson theory. 
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TABLES 



weight 


X = 0.04 


X = 0.07 


X = 0.1 


as 


1 


1 


1 




-1.8 X 10"^^ 


-2.5 X IQ-^^ 


-5.3 X 10-^5 




-1.4 X 10"^"^ 


-1.8 X lO"^'' 


-3.4 X 10"^^ 


ad 


-1.3 X 10-^5 


-8.9 X 10-15 


-3.2 X 10-15 



TABLE L 
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wcnglit; 


./• = 0.07 


j: = 0.13 


./• = 0.19 




1 


1 


1 




-5.7 X 10-15 


-1.3 X 10-20 


6.1 X 10-21 




-2.4 X 10-1'' 


1.2 X 10-21 


-1.3 X 10-21 




-1..3 X 10^1'' 


-3.5 X 10-27 


-2.9 X 10-21 


0-3,12 


9.2 X 10-20 


1.8 X 10-22 


4.9 X 10-18 




-6.7 X lO-^'^ 


1.2 X 10-^^ 


6.2 X 10-25 




2.2 X 10-27 


6.6 X 10-28 


-5.9 X 10-31 


ad,12 


-4.3 X 10-26 


-2.2 X 10-3^ 


-3.5 X 10-31 


as,2i 


5.3 X 10-22 


1.8 X 10-22 


4.9 X 10-18 


Opa:,2l 


-9.3 X 10-29 


1.2 X 10-^3 


6.1 X 10-25 




3.1 X 10-29 


6.6 X 10-28 


-5.8 X 10-31 


ad,21 


-2.5 X 10-28 


1.4 X 10-29 


1.4 X 10-29 


as,22 


1 


1 


1 


ap^,22 


1.0 X 10-20 


7.2 X 10-21 


-1.1 X 10-20 


apy, 22 


1.2 X 10^20 


-5.2 X 10-22 


-6.2 X 10-21 


Od,22 


7.7 X 10-21 


-2.0 X 10-28 


1.6 X 10-2^^ 



TABLE II. 
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wcnglit; 


./• = 0.07 


j: = 0.13 


./• = 0.19 




1 


1 


1 




-4.0 X 10-^5 


-7.2 X 10-21 


1.1 X 10-20 




-4.4 X 10-^5 


5.2 X 10-22 


6.2 X 10-21 




-1.5 X ;u r "^ ' 


1.9 X IQ--^ 


-l.G X 10^27 


0-3,12 


-4.9 X IQ-i^ 


1.6 X 10-15 


6.8 X 10-22 




-2.7 X 10"^^ 


2.3 X 10^15 


3.2 X 10-2^^ 




-4.9 X 10-1^ 


1.6 X 10-22 


-3.7 X 10-34 


ad,12 


-1.3 X 10"^ 


-8.3 X 10-^ 


-1.4 X 10-33 


as,2i 


-4.7 X 10-^^ 


7.5 X 10-11 


6.8 X 10-22 


Opa:,2l 


2.6 X IQ-i^ 


-1.9 X 10-1° 


3.2 X 10-27 




1.7 X 10-1^ 


2.3 X 10-23 


-3.7 X 10-34 


ad,21 


-1.3 X 10"^ 


-8.3 X 10-^ 


-1.4 X 10-29 


as,22 


-1 


-1 


-1 


ap^,22 


7.1 X 10-^5 


2.5 X 10-1^ 


6.3 X 10-21 


apy, 22 


3.1 X 10"^'' 


1.2 X 10^20 


2.7 X 10^21 


Od,22 


2.5 X 10"^^ 


1.0 X lO-i'' 


-1.9 X 10-2^^ 



TABLE III. 
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eigenvector 


u/t 


Lo/t 






{x = 0.04) 


{x = 0.07) 


[x = 0.1) 


s 


0.0111 


0.0089 


0.0063 


s 


0.7221 


1.1180 


1.5592 


Px, Py 


0.7236 


1.1186 


1.5594 


d 


0.7238 


1.1187 


1.5595 


S 


1.4306 


2.2249 


3.1108 


Px5 Py 


J. .^O-LU 






S 


2.5769 


4.0157 


5.6212 


Px, Py 


2.5773 


4.01585 


5.62126 


d 


2.5775 


4.01594 


5.62129 



TABLE IV. 
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eigenvector 


(x = 0.07) 


Lo/t 
{x = 0.13) 


Lo/t 
{x = 0.19) 


8(^0,1(^/3,1 + Sa,2^/3,2), 
s(5q,1(5/3,1 - Sa,2Sf3,2) - ad{5a,lSf3,2 + 5a,2(^/3,l) 


0.0097 
(a = 1.3 X 10"^) 


0.0080 
(a = 8.3 X 10-6) 


0.0054 
(a < 10-10) 


s{Sa,lS/3,l +Sa,2Sl3,2), 
s((Ja,l<5/3,l - Sa,2S/3,2) - ad('^a,l5/3,2 + Sa,2^l3,l) 


0.3494 

(a = 4 X 10-3) 


0.5899 

(a = 1.4 X 10-3) 


0.8935 
(a = 4.4 X lO-'') 


Px(<5a,l5/3,1 + 5a,2<^/3,2), 

Py(<^a,l<^/3,1 + <^a,2^A2), 
Px{Sa,lS/3,l - Sa,2S/3,2) + aPx{Sa,lS/3,2 + 5a,2'5/3,l), 
Py(<5Q;,l<5/3,l - Sa,2Sp,2) - aPy{Sa,lSf3,2 + 5a;,2<5/3,l) 


0.3510 
(a = 3.6 X 10-2) 


0.5905 
(a = 1.3 X 10-2) 


0.8936 
(a = 4.1 X 10-3) 


d('^Q,l<5/3,l + <5a,25/3,2), 
d(5a,l(^/3,l - ^a,2^/3,2) " as(5Q,,i(5/3,2 + ^a,25/3,l) 


0.3512 

(a = 3.4 X 10-1) 


0.5906 

(a = 1.2 X 10-1) 


0.8997 

(a = 3.9 X 10-2) 


s{Sa,lS/3,l +Sa,2Sl3,2), 

s('5„.l'^^,:i,l - <lx.2(\i:2) - «d('\i,l(^.i,2 + '^(i.2<^~,j,l) 


0.6910 

(a = 2.5 X IQ-i) 


1.1707 

(0 = 1.8 X 10-') 


1.7809 

(0 < 10-6) 


Px('^a, 1(^/3,1 + ^a,2Sf3,2), 

Py('^a,l<5/3,1 + (^0,2^/3,2), 
Px(5a,l'^/3,1 - ^0,2(^/3,2) + aPx((^a, 1(^/3,2 + '^a,2(^/3,l), 
Py(<^a,l^Al ~ <^a,2^^,2) - oPy ('^a,l<^/3,2 + (^0,2^/3,1) 


0.6920 

(a = 2.6 X 10-^) 


1.1711 

(a = 2.0 X 10-5) 


1.7810 

(a < 10-6) 


S((5a,l(5/3,1 +<5a,2 (5/3,2), 
s((Ja,l(5/3,l - (5a,25;3,2) " ad{Sa,lSi3^2 + ^a,2(5/3,l) 


1.2030 

(a = 9.7 X 10-2) 


2.1015 

(a = 3.4 X 10-2) 


3.2154 
(a = 1.1 X 10-2) 


Px((5q,1(5/3,1 + ^a,2^l3,2), 

Py((5a,l<5/3,1 +(^0,2^^,2), 
Px(()a,l<)/3,1 - ()a,2<5/3,2j + (jPx(()a,l<)/3,2 + <)a,2()/3,l j, 
Py((5a,l(5^,l - <5a,2(5/3,2) - oPy ((5a,l<5/3,2 + Sa,2Sp,l) 


1.1203 
(a = 1.7 X 10-1) 


2.1017 
(a = 6.4 X 10-2) 


3.2155 
(a = 2.1 X 10-2) 


d((5a,l(5/3,l + 60,2^/3,2), 
d(5a,l(5/3,l - (5a,2(5/3,2) - Os(5a,i(5/3,2 + (5a,2(5/3,l) 


1.1236 

(a = 3.4 X 10-1) 


2.1018 

(a = 1.2 X 10-1) 


3.2155 

(a = 3.9 X 10-2) 



TABLE V. 
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FIGURE CAPTIONS 

Fig. 1 The schematic spectrum of the holon pairing state at temperature (a) above 
and (b) below Tc- 

Fig. 2 The lowest order t-matrix starting from an 1-th orbital state Wi{k) to an interme- 
diate state of momenta ki and — ki. 

Fig. 3 The treated t-matrix with the inclusion of exchange contribution. 
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E 

== s,p,d 

— s,p 
== s,p,d 

(a) T > Tc 



FIGURES 

E 

== s,p,d 

— s,p 
== s,p,d 

s 

(b) T < Tc 

FIG. 1. 
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